Effective range from tetramer dissociation data for cesium atoms 
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The shifts in the four-body recombination peaks, due to an effective range correction to the zero- 
range model close to the unitary limit, are obtained and used to extract the corresponding effective 
range of a given atomic system. The approach is applied to an ultracold gas of cesium atoms 
close to broad Feshbach resonances, where deviations of experimental values from universal model 
predictions are associated to effective range corrections. The effective range correction is extracted, 
with a weighted average given by S.OiO.SRyj-ryy, where R VC [W i s the van der Waals length scale; 
which is consistent with the van der Waals potential tail for the Cs2 system. The method can 
be generally applied to other cold atom experimental setups to determine the contribution of the 
effective range to the tetramer dissociation position. 
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I. INTRODUCTION 

The fundamental quantum properties of atomic sys- 
tems are under intense investigation in recent cold-atom 
experiments, in which long-time predicted universal as- 
pects of few-body physics can be probed. In this re- 
gard, the increasing number of three-body bound states 
(trimers), emerging when the absolute value of the two- 
body scattering length is moved to infinity, known as 
Efimov effect [JJ, was confirmed in different cold-atom 
laboratories [Ufa]- 

Actually, by extending to more than three particles the 
search for fundamental quantum behavior, in the realm 
of four-body physics it is of interest to test the positions 
of four-atom resonant recombination peaks [6-9J, which 
appear as the two-body scattering length is tuned close 
to an s— wave Feshbach resonance (see e.g. [TU]V 

On the theoretical side, the model independence of 
tetramer properties, when few physical scales are fixed, 
was recently numerically established by our group [Til 
[12] , through the analysis of observed correlations between 
tetramer observables, by considering a renormalized zero- 
range (RZR) approach. These correlations were realized 
through appropriate scaling functions, from where one 
can verify that the RZR results are consistent with sev- 
eral other model calculations [TBHTB] . 

Despite the success of the model independence of the 
tetramer scaling functions, the experimental data for ce- 
sium atoms [51 [8, ,9; are shifted with respect to the corre- 
lation function relating the peaks of successive tetramer 
resonances for a zero-range force [TTJ. Among the pos- 
sible effects to explain such a deviation from the results 
obtained with zero-range two-body interaction, it is nat- 
ural to think in a way to introduce the correction due 
to a finite effective range in the model, as suggested in 
Ref. PI]. 

The experimental data for three-body cold-atom re- 
actions show that Efimov resonances and recombination 



minima survive in the intermediate regime when the 
scattering length a exceeds the characteristic interaction 
range by a relatively small factor [T5] . But the positions 
are shifted away from the naive predictions based on the 
universal theory [14[ 115] and finite range potentials have 
been used [111 HO] to address such issues. On the other 
hand, the quantitative reproduction of the experimental 
findings for the positions of the peaks of the tetramers 
resonances has been not yet addressed in detail. 

Within our aim in clarifying the shift in the experi- 
mental data with respect to the zero-range predictions 
for the positions of the tetramer resonances at the disso- 
ciation point, we present an approach on how to intro- 
duce the effective range correction in lowest order in the 
zero-range four-body theory. We obtain the correction to 
the correlation of the positions of two successive resonant 
four-boson recombination peaks. The associated scaling 
function is not fully determined by the trimer properties 
and requires a four-boson scale |llj . Furthermore, we are 
showing by our calculations that the effective range cor- 
rection is not completely parameterized by the trimer and 
tetramer short-range scales, as one could naively expect. 
The contribution of a non-zero effective range affects 
the observed scaling plots, by moving the correspond- 
ing positions in relation to the zero-ranged ones (see e.g. 
Refs. [19l [21] in the case of trimers). Here, it will be 
shown quantitatively how the effective range affects the 
scaling between negative scattering length ratios, corre- 
sponding to the dissociation positions of a trimer and 
related tetramers (for two nearby tetramer states, asso- 
ciated with a given trimer state). 

Relevant to the understanding of the effective range 
contribution to the resonance positions, we should recall 
the nice and clear discussion by Chin et al. [TO] on the 
role of broad (entrance-channel dominated) and narrow 
(closed-channel dominated) Feshbach resonances. As dis- 
cussed, a Feshbach resonance strongly dominated by the 
entrance channel allows a description of a two-atom sys- 
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tern in terms of a single-channel short-ranged model with 
a van der Waals tail. The broad s— wave Feshbach res- 
onance of cesium [HI El E] , where tetramer dissociations 
at the four-body continuum threshold were observed, is 
dominated by an entrance-channel. In this case, when 
the strength parameter is very large (s res » 1), the 
single-channel description is applicable [10] and conse- 
quently, the effective range should be found around the 
values obtained for van der Waals like potentials, tuned 
to the scattering length of the resonance. Indeed, as we 
are going to detail in the present work, the estimated val- 
ues we have obtained for the effective ranges, from the 
positions of tetramer resonances, are somewhat close to 
the corresponding effective range of that class of poten- 
tials. Within our approach to extract the effective range, 
we have considered broad Feshbach resonances appear- 
ing in ultracold gas of cesium atoms, leading to positive 
effective ranges. In contrast, narrow resonances, as dis- 
cussed by Petrov [55], should lead to negative effective 
ranges. In the case of systems with narrow resonances 
and denser spectrum, we should also mention the recent 
work by S0rensen et al. [2"3"] . 

The present work is organized as follows. In Sect. 
|ri] the formalism for two, three and four-boson systems 
is presented, where we adressed the effective range ex- 
pansion of the two-body scattering amplitude, the sub- 
tracted Skorniakov and Ter-Martirosian (STM) equation 



and Faddeev-Yakubovsky (FY) equations. In Sect. [Ill 
the shifts in the position of four-body resonant loss peaks 
from the solution of the STM and FY equations with the 
effective range correction are given. In Sect 



IV 



we per- 
formed an analysis of the positions of the four-boson res- 
onance peaks for the cesium data close to wide Feshbach 
resonances. The effective range, shown to be compati- 
ble with experiments, are compared to the one obtained 
from van der Waals like potentials for the Cs2 system. 
In Sect. [Vlwe present our conclusions and perspectives. 



II. FORMALISM FOR TWO, THREE AND 
FOUR-BOSON SYSTEMS 

A. Two-body scattering amplitude with effective 
range correction 

The s— wave two-boson scattering amplitude input 
to the subtracted forms of the trimer Skorniakov and 
Ter-Martirosian equation [24] and the tetramer coupled 
Faddeev-Yakubovsky equations [25] for negative energies 
e, with the on-shell momentum given by k = \y/—e, can 
be written as 



r(e) 



1/(2tt 2 ) 



ro 
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(1) 



in the lowest order of the effective range ro. In the above, 
the pole of the scattering amplitude is fixed by the bound 



inverse of the scattering length, up to the same order, is 



(2) 



The two-body amplitude ([!]) provides a correction to the 
STM and FY equations in leading order in ro, and as we 
verified the subtracted forms of these equations are able 
to provide finite results dealing with the extra power of 
momentum in the amplitude. 

In the following, we present some details on the formal- 
ism for the subtracted equations we have used, together 
with our approach to find the effective range correction 
for the position of the tetramer resonances. 



B. Subtracted STM equation 

The Faddeev components for three-boson bound states 
for the zero-range potential can be rewritten in terms of 
spectator functions, which are solutions of a subtracted 
Skorniakov and Ter-Martirosian (STM) integral equa- 
tion, given bv|24j: 



\K 



ij,k . 



K 



ik,j , 



(3) 



where ey ^ is the (ij) subsystem energy in the three-body 
system. Note that we are going to introduce the two- 
boson amplitude carrying the effective range correc- 
tion in Q. The subtraction of the free 3-body resolvent 
regulates STM equation, and projected operator, is: 
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(+) or virtual state (— ) energy, k 



e 2 l The 



where Hq is the three-body free hamiltonian and — /x§ is 
the subtraction energy scale. The form factor appearing 
in the two-body T- matrix of the pair (ij), in the relative 
momentum id given by (pylxij) = 1- The three- 
body regularization parameter fi^ keeps under control the 
Thomas-collapse, and determine the three-boson observ- 
ables, and can be parameterized by one s— wave physical 
quantity (see e.g. [H]). The subtracted STM equations 
was also used to calculate the position of triatomic con- 
tinuum resonances for large negative scattering lengths 
[37] . which were observed by the Innsbruck group [2J, 
through three-atom recombination peak in a cold cesium 
gas close to a Feshbach resonance. 



C. Subtracted FY equations 

The tetramer energies are found by solving the 
Faddeev-Yakubovsky (FY) equations properly regulated 
in the limit of the zero-range interaction, which are writ- 
ten as a set of coupled subtracted integral equations 
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[331 [HHS] given by: 



extended to include the effective range, as 



l£«,fe) = 2r ( e y,fc) 
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The subtractions in the kernel come at the level of the free 
4-body propagators in Eqs. ([5| and ([6]), which are used to 
regulate the FY equations for the contact potential. The 
two-boson scattering amplitude (J7J), which are introduced 
in the FY equations above, carries the effective range 
correction in the calculations that are presented in the 
next section. 

The projected 4-body free resolvent, are sub- 

tracted at an energy — //|: 



S$ ik = (Xij\[[E-H { 
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where H is now the free four-body Hamiltonian. Notice 
that when introducing Q^\^, eq. (jij, in the FY equa- 
tions, the free three-body hamiltonian should be substi- 
tuted by the four-body one. The energy of the two-body 
subsystem (ij) appearing as arguments of the two-boson 

and 
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ij,k 



scattering amplitude in Eq.s (|5j) and 
£ij,kh associated with a virtual pair in trie 3+1 partition 
and in the 2+2 partition, respectively. 

The subtraction scale — //§ in Q^\ k the right-hand side 
of Eq. ([5]), fix the trimer properties consistently with the 
STM equation, but independently on the scale /X4. Look- 
ing closer to this term and taking into account the left 
hand side, the position of the trimer pole is guarantee to 
be same as the bound state energy obtained by solving 
the subtracted STM equation. The other terms in the 
FY equations, in ([5| and ([6]), are regularized with an in- 
dependent subtraction scale — This method ensures 
that the three-body subtraction scale fixes the three- 
body properties consistently with the subtracted STM 
equation, whereas the regularization parameter \x A fixes 
the four-body observables [2"5] . 



III. POSITION OF FOUR-BODY RESONANT 
LOSS PEAKS 



The four-atom recombination resonates when the 
tetramer has zero energy and placed at the four-body 
continuum threshold. In cold atom traps the resonance 
is reached by tuning the negative scattering length close 
to a Feshbach resonance (see e.g. [H]). The theoretical 
description of the scaling function, which correlates the 
values of the negative scattering lengths where successive 
tetramers reach the continuum, introduced in is now 




(5) where a N is the position of the peak of the three-atom 
resonant recombination for a < 0, and N is the scat- 
tering length for which the excited iV-th tetramer dis- 
sociates and meets the four-body continuum. The solu- 
tions of the subtracted STM Q and FY equations fl5pj ) 
depend only on a, ro, and on the three- and four-body 
scale parameters. Then for the zero energy tetramer, the 
value of the negative scattering length depends only on 
ro and the short-range three- and four-body momentum 
scales at the subtraction points. The three-body scale 
can be parameterized by ijv 3 , while the four-body one 
by a^ 3 N in the determination of ajf 3 N+1 - In this way 
^ can be built. Furthermore, the scaling function A is 
calculated with the two-boson scattering amplitude ([!]) 
and the solutions of the FY equations, as well as of the 
STM equation, expanded up to first order in ro. 

The dimensionless function A is determined by solving 
the subtracted STM [Mj and FY [55] equations consid- 
ering the expansion of the atom-atom s— wave amplitude 
up to order k 2 as given by ([I]) and the corresponding 
relation between the scattering length and virtual state 
energy The solutions of the subtracted FY equations 
depend only on a, ro, and on the three- and four-body 
scale parameters. Then, for the zero energy tetramer, the 
value of the negative scattering length depends only on 
r and the short-range three- and four-body momentum 
scales at the subtraction points. The three-body scale 
can be parametrized by , when solving the STM equa- 
tion with ro included, while the four-body scale is corre- 



lated to 



N 3 ,N 



in the determination of 



N 3 ,N+1' 



Then, 



the scaling function <[8j) is built without reference to the 
scale ratio, by eliminating this dependence in the ratio 



l N 3 ,2 



/a N in terms of ajj 3 1 /a 



In our calculations we have only obtained results for 
tetramers below the ground state trimer, which should 
be enough for our study, as the values of the scattering 
length are much larger than the short-range length scales, 
corresponding to the three- and four-body subtraction 
points. The scaling function A (a[ i/a^, ro/a^), which 
provides the correction to the position of the four-atom 
resonance is shown in the 3D plot of Fig. [TJ We have 
presented results up to r /a^ ~ 0.5, for 0.40< a\ l /a[ < 
0.51, considering the relevant region where recombina- 
tion data for cesium have been measured. In this case, 
and for the validity of the expansion in the 



ro 



n,ij 



effective range one should have r << 1 , and the ex- 
pansion of the scattering length ^ is questionable. For- 
tunately, we found that the coefficients for the effective 
range correction are fairly small, allowing to extend the 
validity of the expansion to the region where data are 
found. 

Within the region of the plot given in Fig.[TJ the surface 
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FIG. 1. (color online) The scaling function A expressing the 
dependence of aj 2 / a i m a i,i/ a i an d ro/a^. Calculations 
are performed in a leading order expansion of the solutions of 
the FY equations with the two-body amplitude see text 
for further explanation. 



can be parameterized by 



A(x, y) = c ™» ( x " 0A5 ) m y r ' 



(9) 



where x = a{ l /a l , y = ro/\a 1 |, and (m,n) > 0. The 
coefficients are given in the following Table |TJ 

TABLE I. Coefficients for the parametrization |9|. 

COO Cio Coi Cn C20 CQ2 



0.932 0.724 -0.144 0.347 -0.645 0.001 



The negative value of coi decreases the function A for 
positive effective ranges, which suggests that the region 
where the two-body amplitude increases with respect 
to the zero-range one, for a < and positive effective 
ranges, is important for the tetramers. In this case, larger 
variations of the scattering length toward negative values, 
compared to the trimer value of aj~ , are required to tune 
the tetramers to zero energy. 

Two theoretical results are used for comparison with 
the results of the scaling function A, using the plot of 
Fig. the first one by Deltuva [IS] (0.4435,0.8841), 



which quotes tq/ 



0.33 



compared to 0.29 



extracted from the plot, and another one by von 
Stecher Q3], (0.45,0.88), with r /\a^\ =0.38 [27] com- 
pared to 0.36 obtained from the plot. This compari- 
son with completely different models, separable poten- 
tial model p~6] and local Gaussian potential [13] , gives us 
confidence on the universality and utility of the function 
A, to analyze the data. This comparison suggests a lower 
bound of |ro/a^~| > 0.02 in our extraction method. 
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Although the expansion (J9]) contains nonlinear terms 
in r , it is built in linear order in the solutions of the 
subtracted STM and FY equations. The second order 
terms are induced by the curvature of A{x, y) in x and y. 
In our case the calculation is performed for the ground 
state trimer and two successive tetramer states. 

The coefficients for the effective range correction are 
expected to be smaller than the ones associated with 
the variation of aJi/a^. Indeed, the expansion around 
a i i/ a i ~ 0A5 shows that the linear coefficient for ro/aj" 
correction gives coi/coo = —0.15, while cio/coo = 0.78 is 
5 times larger. The amazing smallness of coi and C02 
with respect to the other coefficients, reflects that a frac- 
tion of the effect from the effective range is absorbed by 
the variation of the short-range four-body scale, but not 
all. Given that, it is quite obvious that the correction 
coefficients for the effective range contribution, shown in 
Table [I] should be somewhat smaller than the coefficients 
associated with the expansion parameter a^ 1 /a^ . 



FIG. 2. (color online) Lower-order effective range correc- 
tion for the scaling plot a\ 2 / a \ = A(a[ 1 /a^ , ro/a^) for 
the correlation between the positions of the peaks for succes- 
sive tetramer resonances at the continuum threshold, com- 
pared to other theoretical results and experimental data. We 
present our results for different ro intervals (the upper solid 
line is from jTT) compared to a few recent calculations as 
indicated inside the frame [Refs. [13] (triangles), [16] (dia- 
monds) and [28] (square)]. The experimental data of Refs. [6], 
[9] and [8] are also indicated inside the figure. The lowest 
band (gray) corresponds to ro for the van der Waals like 
potential for cesium atoms with values within the interval 
0.4 < —ro/a^ < 0.5 for a reference value of ~ — 
(see Table M. 

In Fig. [2| the 2D plot of our results for the surface 
A(x,y), corresponding to the scaling function ([8]), is 
shown. The bands in the plot gives the shift of the scal- 
ing function due to the finite values of r / 1 1 in the 
intervals depicted in the figure. For our reference we also 
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plot the theoretical results from [T3] and [TB], and the 
experimental values from Refs. [SI [5J [9] . The data from 
Ref. 7] have large error bars and suggests r ranging 
from negative to positive values. 



IV. ANALYSIS OF THE CESIUM DATA 

The effective range extracted from the experimental 
data is shown in Table [TTJ They were obtained just by 
inspection of Fig. [2j by comparing the experimental val- 
ues with the scaling curve including the effective range 
correction. The shift of the data yields values of ro/aj - 
larger than the lower bound for the theoretical extrac- 
tion. We didn't present results for [7] as it has larger er- 
rors. We just notice that the effective ranges from these 
data varies from negative to positive values. The ex- 
traction of ro from the data for cesium atoms close to 
Feshbach resonances are within errors consistent. The 
extracted values of r vary from 2 to 5 times the van 
der Waals length (-Rydw)' wrnc h f° r the Cs 2 system is 

r^ w = loi-ooo ma. 

The effective range of potentials with van der Waals 
tail —Cs/r 6 at large distances has the approximate for- 
mula given by [29H5T] : 



-I^"-'- (t)"*(!-0". « 

where a = 2vr/r(l/4) 2 i? vdw ^ 0.955978 i? vdw is the 
average scattering length of van der Waals potentials [31] , 
with the van der Waals length i£y d w ~ ^mCe/h) 1 ^- 
Considering the van der Waals length one gets that a = 
96.5 do, and 

r ~ 2.7894 i? vdW = 2.9179 a <d ~ 281.7i?a . (11) 

In the actual cold cesium experiments listed in Table [TTJ 
the factor d varies from 1.26 (for the second tetramer res- 
onance) to 1.55 (for the first tetramer resonance), which 
gives ro within 3.5 — 4.3i£ vc m7. The corresponding in- 



terval of 0.4 < 



< 0.5 for a reference value of 



a-i ~ —dRy^yr, can be identified with the lowest band 
(gray) shown in Fig. [2] 

The estimated values of r extracted from the shift 
of the data with respect to the zero-range calculations, 
shown in Fig. [2j are given in Table [TTJ The values for 
the effective range are found between 2/?^^ and ~ 
5i? vd ^, while Eq. (Ill applied to the set of cesium 
data suggests the interval 3.5 — 4.3i? vd ^ represented 
by the lowest (gray) band in the figure. The present ex- 
perimental errors overlaps the cesium data, as one verify 
in Fig. [2] By performing the error weighted average of 
the extracted effective range values from the data given 
in Refs. [8] and [9] we obtain 3.9±0.8 Ry^w ^ n particu- 
lar, the data set {0.47(1), 0.87(1)}, provides an effective 
range of about 4.8i? vd ^, which within la is consistent 
with the values of the van der Waals potential. The data 



set {0.47, 0.84} [B] with no quoted errors suggests a large 
value of r Q > 5 -R vd \y- 



TABLE II. Extracted effective ranges from the shift of the 
experimental peaks of the four-atom losses for a gas of cold 
cesium atoms close to a Feshbach resonance, using the scaling 
plot shown in Fig. [2] The errors in the extracted values of ro 
are estimated from the figure as well. 
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It is noticeable that the dissociation point of ce- 
sium ground state Efimov trimers for negative scatter- 
ing lengths are found in a narrow band (8], independently 
on which Feshbach resonance the system is tuned (see 
Table [TTJ) . Theoretical works [231 I52H3"5] addressed the 
interesting issue of the physical mechanism for the dom- 
inance of the two-body properties on the position of the 
first Efimov resonance for a < 0, when the trimer meets 
the continuum. The results in Table [TTJ show a range 
of values for ro extracted from the tetramer and trimer 
resonances, which within errors (excepting one data set) 
are consistent with the effective range values from the 
van der Waals type potential. The numerical analysis of 
Wang et al. [32] has shown that two-body interactions, 
which suppress efficiently the wave function for separa- 
tion distances less than r , have a± ~ — 9i? vd ^r, deter- 
mined by two-body properties, with this class of systems 
closely related to entrance-channel dominated Feshbach 
resonances. The dominance of the single channel poten- 
tial with the van der Waals tail in the three-cesium reac- 
tions, with suppression of the short distance wave func- 
tion, putted forward the narrow band where the trimer 
dissociation position is found. It also suggests that the 
effective range seen through the shift of the position of 
the four-atom recombination peak with respect to the 
zero-range results, should be given by the van der Waals 
tail of the potential, as well. 

It is worthwhile to address the validity of the correc- 
tion due the effective range, considering the main results 
presented in Fig. [2] and Table [TT] In view of the expansion 
given in Eq.Q, which at a first glance would be doubtful, 
when ro/|a^"T~ af 2 /cii , we found consistence between 
the present calculations with the results obtained with 
short-range separable[5B] and local[17] interactions. The 
leading coefficient for the expansion in ro/|a^~ | is one fifth 
of the corresponding coefficient for the dependence of the 
scaling function in a^/aj - (see Table[lj), which shows that 
the extracted values of r using the scaling plot comes 
from small corrections to the zero-range scaling function. 
Therefore, as observed in Fig. [2j the agreement of the 
expansion we have considered can be extended also to 
cases near af 2 /o,i ~ r o/|ar|. 
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V. CONCLUSION AND OUTLOOK 

In this work, we present an approach to extract the 
effective range of a given atomic system from the shifts 
in the four-body recombination peaks, at the tetramer 
dissociation threshold, with respect to the zero-range re- 
sults. We solved the trimer and tetramer subtracted zero- 
range integral equations with two-body amplitudes car- 
rying the effective range in lowest order. Besides, the 
trimer and tetramer subtracting scales, no further pa- 
rameters are introduced in the calculations. The correla- 
tion between the negative scattering lengths where suc- 
cessive tetramers dissociates at the four-body continuum 
is extended to include the correction due to the effec- 
tive range, which is model independent checked against 
different short-range potential model calculations. We 
found that the effective range correction of the on-shell 
two-body amplitude in the calculation of the trimer and 
tetramer dissociation points, when presented as a cor- 
relation plot of the successive positions of the tetramer 
resonances is not completely parameterized by the trimer 
and tetramer short-range scales, as one could naively ex- 
pect. This unexpected property turned to be essential to 
single out the effective range from the shift in the posi- 
tions of the four-atom recombination peaks with respect 
to the zero-range theory. 

We applied our proposal to the cesium data mea- 
sured at broad s— wave Feshbach resonances [5J [H]- 
In our analysis we considered the shifts of the observed 
recombination peaks at the tetramer dissociation posi- 
tions with respect to the zero-range results. The effec- 



tive ranges were found within the interval 2 i? vc j^ ^ 
r o S 5-R v dw w ith a reference value of 3.9±0.8 -R vc jw, 
obtained by the error weighted average of the r 's ex- 
tracted from the data of Refs. [8,9. These values are 
consistent with the van der Waals potential tail for the 
broad s— wave Feshbach resonances of the Cs2 system, 
providing r within the interval of 3.5 — 4.3 R V( ^\^ for 
the scattering lengths where the resonances were found. 
Our results put a strong evidence on the prevalence of the 
entrance-channel dominance in the physics of the univer- 
sal tetramers formed with broad Feshbach resonances. 

The question on how the effective range is sensible 
to Feshbach resonance manipulations in the scattering 
length as suggested by pTj ) has to be considered. Once 
improved experimental data are available, it will be pos- 
sible to verify in more detail how Feshbach resonances 
from different setups affect the ro values extracted from 
the trimer and tetramer dissociation points. The present 
analysis, based on the universality of effective range cor- 
rection to the position of tetramer resonances, can be 
generally applied to any other cold-atom bosonic system 
close to Feshbach resonances. 
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